How many atoms get excited when they decay? 
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We study the dynamics of an atomic two-level system decaying by spontaneous emission of light. 
Subject to continuous detection of the radiated field, the system tends with certainty to the ground 
state in the long time limit, but at initial times the excited state population exhibits non-trivial 
stochastic behavior. Employing methods from Ito calculus, we characterize this behavior, and we 
show, for example, that the emitter, as a result of in-phase homodyne measurements, may become 
fully excited during the decay process while heterodyne and out of phase homodyne measurements 
do not drive the atom completely into the excited state. 



I. INTRODUCTION 

With the ability to control single emitters in the form 
of trapped ions, atoms, NV-centers, Josephson qubits, 
etc., the real time dynamics of single quantum systems 
subject to measurement constitutes a very active research 
field vifith wide perspectives for fundamental investiga- 
tions as well as for quantum metrology, quantum com- 
puting, and quantum feedback and control. 

In quantum optics light emitting systems are often 
studied by subjecting the radiation to photon counting, 
heterodyne or homodyne measurements. In case of con- 
tinuous measurements, the conditioned quantum state 
of the emitter evolves according to a stochastic master 
equation (SME) or, if the state remains pure as a result 
of the measurements, a stochastic Schrodinger equation 
(SSE) [H-lll. The SME and SEE exphcitly reveal the fiuc- 
tuation and noise characteristics of the emission process 
and the detection back-action, while averaging over an 
ensemble of measurement histories leads to the conven- 
tional Lindblad master equation describing the state of 
an unobserved emitter subject to decoherence and decay. 
Dynamical equations similar to the measurement SSEs 
also arise in the description and efficient simulation of 
open systems [1, @ as well as in collapse models aiming 
to shed light on the measurement problem 0] . 

In this paper we study the state of a decaying two- 
level quantum system subject to continuous, unit effi- 
ciency heterodyne or homodyne detection of the emitted 
light. For photon counting it is well known that the state 
of a two-state emitter experiences continuous decay in 
periods with no photodetection events, while a photon 
count event is accompanied by a quantum jump of the 
system into the ground state [s!]. For heterodyne and 
homodyne detection, on the other hand, the continuous 
measurements of the intensity interference signal between 
the source field and a local oscillator field cause the ex- 
citation of the emitter to undergo random, diffusive time 
evolution d, |^. Such behavior is illustrated in Fig. [U 
where the excited state population (solid curve) of a two- 
level emitter under homodyne detection is compared with 
the exponential decay (dotted curve) of the unobserved 
system. In this article, we use techniques from Ito cal- 



culus to analytically characterize different aspects of the 
measurement-induced dynamics of the light emitter. In 
particular, we determine the probability that the emitter 
during its decay attains any specified degree of excita- 
tion, and we calculate the mean time the observed atom 
spends above its initial excitation as well as the aver- 
age time needed for the system to reach different excited 
state populations below the initial value. 




FIG. 1. A single realization of the upper level population 
Ct (solid curve), exhibiting transient fluctuations that exceed 
the initial value Co = 0.5 during homodyne detection. The 
dotted curve shows the exponentially decaying, unconditioned 
excited state probability pee{t) with pee(O) — Co- 

In Sec. ini we introduce the Lindblad master equation 
and the stochastic wave function descriptions associated 
with photon counting and heterodyne and homodyne de- 
tection. For heterodyne and homodyne detection, we de- 
termine the resulting stochastic differential equations for 
the excited state population of the emitter. In Sec. IIIIl 
we present numerical simulations of homodyne detec- 
tion with arbitrary local oscillator phases. Focusing in 
Sees. IIVI and |V] on heterodyne detection and homodyne 
detection with the local oscillator initially in phase with 
the emitter dipole, we derive adjoint Fokker-Planck-type 
equations for interesting probabilistic and average prop- 
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erties of the conditioned atom. In Sec. IIVI we thus de- 
termine the probability that the excited state population 
reaches the upper endpoint of an interval before the lower 
one, and we heuristically discuss the attainability of zero 
and unit excitation. This permits analytical evaluation of 
the probability for the excited state population to reach 
any level. In Sec. |Vl we obtain an expression for the av- 
erage time it takes the upper level population to leave 
definite intervals, and we use this result to rigorously 
discuss the attainability of zero and unit excited state 
occupation. We also determine the mean time the emit- 
ter subject to detection spends above the initial excited 
state population as well as the average time it takes the 
upper level population to attain a specified value below 
the initial excitation. Sec. IVII concludes the paper. 

II. STOCHASTIC DESCRIPTION OF A 
TWO-LEVEL QUANTUM EMITTER 

We consider a two-level quantum system with excited 
state |e) decaying into the ground state \g) by sponta- 
neous emission of electromagnetic radiation. The time 
evolution is treated in the electric dipole, rotating wave, 
and Born-Markov approximations, and if no measure- 
ments are performed on the emitted field, the system will, 
as a result of the interaction-induced entanglement with 
the quantized radiation field, evolve into a mixed state 
obeying the interaction picture master equation (we set 
h = 1 throughout) [H, i] 

dp{t) - -f{a^pit)a+ - ^{<j+a^,pit)})dt, (1) 

where 7 is the decay rate, (t_ = \g){e\ is the lowering 
operator, and (T+ — \e){g\ = ul is the raising operator of 
the two-level emitter. 

The solution to the master equation Eq. ([1]) yields the 
well known exponential decay of the excited level popu- 
lation and the coherences, 

PeeW = l-PssW -Pee(0)e-T* (2) 

When the light emitted from the decaying system is 
subject to continuous photon counting, homodyne or het- 
erodyne detection, the atom-field entanglement mediates 
a detection back-action on the system state which, in the 
case of unit efficiency measurements, restores the pure 
state character of the emitter, provided it is initially pre- 
pared in a pure state. 

For continuous photon counting Q , the evolution is 
split up into no-jump and jump dynamics, such that the 
pure state \i^{t)) = Ce{t)\e)+Cg{t)\g) of the emitter during 
intervals with no photon detection experiences a damping 
rate of 7/2, and the upper level population Ct = |ce(i)P 
evolves as Ct = Coe~'''*/(l — Cq + Cqc"'^*). Associated 
with the detection of a photon, this non-exponential de- 
cay of the excitation is with probability ^Ctdt interrupted 



at time t by a jump of the state vector into the ground 
state, where the system remains in the absence of an 
external driving field. The average excited state popula- 
tion over many realizations of this dynamics is given by 
Eq. ([2]): jumps occur with the rate jpeeit), and on aver- 
age, a fraction pee(O) of the emitters will reach the ground 
state by a jump, while the remaining fraction Pgg{0) en- 
ter the ground state by continuous, exponential no-jump 
evolution. For both the jump and no-jump dynamics, 
the excited state population of the two-level emitter is 
always lower than in the initial state. 

For homodyne detection with local oscillator phase (/>, 
the system evolves according to the interaction picture 
stochastic Schrodinger equation (SSE) P, 2, 8] 

d\^P{t)) = i7(-cr+(7_ + {cr+e"^ + a^e-'^)a-e~"'' 

+^(a-e-'^ - + ^-e"^*)) m))dWt, 

(3) 

where (. . .) denotes the quantum mechanical expecta- 
tion value {%l){t)\ . . .\%lj{t)) . The independent, Gaus- 
sian distributed Wiener increments {dWt\t>G with mean 
'Ej{dWt) = and variance Var(dWt) = dt represent the 
random outcomes of the homodyne detection through 
the appropriately scaled differential measurement current 
dqt = 7(cr+e*'^ + a_e-''f')dt + ^dWt and obey the Ito 
calculus rules dW"^ = and dWtdt = [2]. 

The corresponding heterodyne SSE reads P, Q 

d\m)^i{-\'y+'^- + (^+>^- - \{<^+){'y-))m))dt 

+ - {<J^))m))dZu (4) 

where the complex Wiener noise dZt = -^{dWt,x + 
idWt^y) is defined in terms of two real-valued, in- 
dependent Wiener increments dWt^x and dWt^y with 
dWt^xdWt^y = 0. dZt thus obeys dZ^ = and dZ^dZt = 
dt. The complex differential measurement current is 
dqt = ^{a-)dt + y^dZt. 

Upon calculating d{\ip{t)){il:{t)\) = d\ip{t)){ip{t)\ + 
\'ip{t))d{'tp{t)\+d\i;{t))d{i:{t)\ fromEq. © or Eq. and 
using the Ito calculus rules, we verify that the ensem- 
ble average F,{\4>{t)){ip{t)\) fulfills the master equation 
Eq. (P). 

Applying the expansion |V''(i)) = Ce(i)|e) -I- Cg{t)\g) 
in Eq. ^ and Eq. (jH), we obtain coupled equations 
for Ce(t) and Cg{t). Using the Ito calculus rules, and 
the normalization |ce(i)p -f- |cg(i)p — 1, we find a 
stochastic equation for the excited state population Ct = 
|ce(t)p. Similarly, by taking the differential of (pt = 
Im[ln(c*(i)ce(i))] , we find an equation for the phase dif- 
ference Pt — feit) — fgit) between Ce{t) — 106(^)16"^"=^*^ 
and Cg{t) = |cg(t)|e''^«(*). 
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For heterodyne detection, the resulting equations read 
dCt = -jCtdt - ■^y/^^Cf{l~Ct)dWt 

= A{Ct)dt + BiCt)dWt (5) 



dipt 



l-Ct 



{iiin{ipt)dWt.^ + coii{ipt)dWt.y) 



where we have defined dWt — cos{ift)dWt^x — 
su\{Lpt)dWt,y, which is seen to be a standard Wiener in- 
crement upon using the properties of dWx,t and dWi 



v.t 



and their independence of ipt ■ Consequently, Eq. ((S]) con- 
stitutes a closed equation on standard Ito form for Ct- 

For homodyne detection we get instead the two cou- 
pled equations 

dCt - -7Ctdt - 2 V7 cos((^i ) ^Cfil~Ct)dWt (6) 



2C? - Ct 



d'ft = 1 ^^^ _ s\TL{2(pt)dt ~ ^s\w{Cpt)\j ^ _ fj^ dWt, 

(7) 



Ct 



where ipt = ipt — <j) and where we observe from Eq. ([6]) 
that the initial phase (po may be restricted to [0, -^J (by 
absorbing a possible minus sign in the definition of dWt). 

For arbitrary local oscillator phases, the time evolution 
of Ct under homodyne detection is generally coupled to 
that of the phase variable (pt- However, if initially the 
emitter dipole is in phase with the local oscillator and 
(po = 0, Eq. ^ shows that = and thus cos((^t) — 1 
for all future times < > 0. Since the magnitude of the 
measurement back-action induced by the diffusion term 
in Eq. ([6]) is largest when | cos{(pt)\ = 1, we will refer to 
this case as optimal homodyne detection and we get a 
closed equation for Ct 



dCt = -iCtdt - 2^^Cf{l-Ct)dWt 
= A{Ct)dt + B(Ct)dWt. 



(8) 



Fig. [T] shows a single realization of Eq. (jS]) found by nu- 
merical simulation using the Milstein scheme 10]. Unlike 
the exponentially decaying excitation of the unobserved 
atom, the excited state population conditioned on the 
optimal homodyne detection signal reaches a maximum 
value well above the initial excitation. It is the purpose 
of this paper to characterize the stochastic dynamics of 
the excited state population. In particular, we shall de- 
termine the probability PT„<oo(y) that the excitation Ct 
of the emitter (in finite time T„ < oo) reaches different 
population levels u above the initial value Cq = y for the 
different detection schemes. 



dipole, the coupled equations ^ and ([7]) for the popula- 
tion and phase variables Ct and (pt cannot be approached 
analytically and in this section we estimate PT„<oo(y) by 
numerical simulations of Eq. ([3]) . 

In Fig. [21 we present data for simulations, where 
we begin with excited state population Co = y = 
0.3,0.5,0.7,0.9, and where the homodyne local oscil- 
lator phase initially differs from the dipole phase by 
(^o/f — 0.1,0.3,0.5,0.7,1. For each initial population 
and phase, we propagate 5000 independent solutions of 
Eq. ([3]) and register the fraction of trajectories Ct that 
exceed different values u above the initial value Co = y. 
This fraction is thus a numerical estimate for the desired 
probability Ptu<oo(j/) that the emitter starting with an 
excited state population y reaches the excitation level 
u > y in finite time. 

For the simulations, we use an explicit 1.0 strong or- 
der Ito- Taylor scheme fl^ with a global pathwise er- 
ror of 0{At), where At is the temporal step size. All 
simulations have been carried out on the time interval 
[07~^, 57""'^] and with a step size of At = 10~*7~^. In 
the case of optimal homodyne detection with ipo = 0, we 
have for all values of u G [0, 1[ verified that our simula- 
tions of Pt„<oo(2/) accurately to within ±0.02 reproduce 
our analytical derivation, Eq. ([T5)) . 

Each simulated emitter achieves with certainty the ini- 
tial excited state population and, depending on the ho- 
modyne detection phase, it reaches with different prob- 
abilities excitation levels above the initial value. The 
results of our simulations are reported in Fig. [5] as sets 
of curves decreasing from unit probability beyond differ- 
ent upper state population levels u. The upper dashed 
curves are exact results for in-phase homodyne detection. 
We see from Eq. ([7]) that an initially vanishing phase 
difference (j)t remains zero for all times, and that this 
choice of phase leads to the strongest stochastic back- 
action in Eq. (jG]). In Fig. [2l this is reflected by higher 
excitation levels during the atomic decay process. For 
(po — f , we see from Fig. [2] that Ct never exceeds its 
initial value during the homodyne probing of the emit- 
ted radiation. For smaller values of the detection phase, 
</5o/f = 0.1,0.3,0.5,0.7, the atom may be become ex- 
cited, but the level of excitation decreases with increas- 
ing phase, and none of our simulations lead to a value of 
Ct exceeding 0.999. 



IV. PROBABILISTIC DESCRIPTION OF 
EXCITATION DYNAMICS UNDER 
HETERODYNE AND OPTIMAL HOMODYNE 
DETECTION 



HI. NUMERICAL SIMULATION OF GENERAL 
HOMODYNE DETECTION 

For general homodyne detection, with a local oscilla- 
tor phase which may differ from the initial phase of the 



In case of heterodyne and optimal homodyne detec- 
tion, the excited state population Ct evolves indepen- 
dently of the phase variable according to the closed equa- 
tions ([5]) and ([5]), respectively. These population equa- 
tions, however, cannot be reduced to linear equations and 
are thus not explicitly solvable by the reduction method 
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FIG. 2. The probabilities PTu<oo(y) that the upper level 
population Ct attains the value u above its initial value 
Co = y — 0.3, 0.5, 0.7 are shown as groups of curves from left 
to right in larger figure while the results for Co = y = 0.9 are 
shown in the smaller figure. Each group of curves presents cal- 
culations carried out with initial homodyne detection phases 
<^o/f = 1,0.7,0.5,0.3,0.1 (from left to right), converging to- 
wards the analytical result for optimal homodyne detection 
with ifo = (dashed curve). 

M- 

Fortunately, Eqs. ([S]) and ^ may be analytically in- 
vestigated using standard results from Ito calculus. To 
begin the analysis, we observe that once Ct hits 0, the 
excitation is irretrievably lost since dCt\ct=o — and the 
system remains in the ground state. If, on the other hand, 
Ct reaches unity, dCt\c t=i = ~ldt, and the atom with 
certainty decays away from the excited state again. With 
these boundary conditions a standard analysis [ill, [13 
implies that Eqs. ([5]) and ([8|) have pathwise unique, (al- 
most surely) continuous, and Markovian solutions on 
[0,1], as also required by the conservation of probabil- 
ity and the Markov approximation used in the derivation 
of the SSEs. 

In order to quantitatively describe the dynamics of 
Ct, we introduce two characteristic times, visually repre- 
sented in Fig. [3l We define the first hitting time T„ (y) of 
u to be the first time that Ct with Co = y G [0, 1] attains 
a given value u in [0, 1], i.e. Tu{y) = inf4>o{Ct — u\ Cq = 
y}. Similarly, we define the first exit time Ta.b{y) o/]a, h[ 
to be the first time that Ct with Cq = y ^ [a, b] reaches 
either a or b, that is Ta^iv) = inft>o{Ct (f. ]a, b[\ Cq = y}. 
We note that Ta^iv) Tain{Ta{y) ,Tb{y)} . 



A. Probability to hit b before a 

To quantify the extent to which the emitter becomes 
excited as a consequence of the measurements, it is con- 
venient to introduce the probability PTb<Taiy) that Ct, 



Ct 



1- 




0-^ ' ' »'t 

Tb{y) = Ta^biy) Ta{y) 

FIG. 3. First hitting times Ta(y) and Tb(y) of the levels a and 
b, respectively, and the first exit time Ta.b(y) of the interval 
]a, b[ for a particular realization of the upper level population 
Ct given in Eq. (0) or Eq. ^ with Co ^ y £ [a, b]. Since Ct 
hits b before a in this case, T(,(j/) < Ta(y) and Ta.biy) ~ Tb(y). 

commencing from Cq = y € [a,b], hits the value b before 
a. Following J^], we shall heuristically obtain a differen- 
tial equation for Px^^Taiy)- 

Let P{y) = PTi<Ta{y), and note that the probability 
that Ct, starting from Cq — y ]a, &[ at time t — 0, 
reaches a or 6 already after an infinitesimal time dt is 
vanishing. Thus, conditioned on the value of Cdt, the 
probability to reach 6 before a at time dt is P{Cdt) and 
we conclude that the average value E(P(Cdt)) over all 
possible realizations of C^t starting at Cq = y must equal 
the initial probability P{y). At the same time, we have 
from Eq. or Eq. ^ and the Ito rules that E(dCo) = 
A{y)dt and E(dC^) = B{y fdt. 

Combining these two results and using that dCo = 
Cdt -Co^Cdt-y,^e obtain 

Piy) = E(P(C,t)) 

= E(P(rfCo+y)) 

= E(p(y) + P'{y)dCo + lp"{y)dC^) 

= Piy) + P'{y)A{y)dt + ]^P" {y)B{yfdt, (9) 

where we have Taylor expanded P{dCo + y) around y 
to second order in dC^ and used the Ito rule dW^ — dt. 
Now, subtracting P(y) on both sides of Eq. ([9]) and divid- 
ing through by dt, we obtain the sought-after differential 
equation for PT^KT^iy) 

LyPn<TAy) = o (10) 

subject to the natural boundary conditions Pt^ktAo) = 
and PTb<Ta (b) = ^, where we have defined 

^y-My)^^ + lBiyr£,, (11) 

For more details see Ch. 15.3 of [l^ and Thm. 4 in 
Ch. 3.15 of [13. 
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The solution of Eq. (|10p can be found explicitly by 
direct integration to be 



Pn<TAy) 

with y G [a, b], where 



S{a,b) 



(12) 



S{u,v)^ I Q{z)dz, Q(z) = e '^'^"^ (13) 



We note that for Q and thus 5 to be a priori well-posed 
we must require that the diffusion term B{Ct) in Eqs. ([5]) 
and © obeys B{Ctf > 0, which is fulfilled for < a < 
b < 1. 



1. Heterodyne detection. 

For heterodyne detection, the integrals in Eq. (|T3| can 
be evaluated analytically, and we obtain 



Q{z) 

S{u, v) 
with limiting behavior 



1 - z 



e \ z / 



S{0,b)<oo, 5(a, l) = oo for a, & G [0, 1[, (14) 

where Ei{z) — J°° ^—j-dz' is the i^i-function with 
Ei{Q) ^ oo and Ei(oo) = 0. 



2. Optimal homodyne detection. 



For homodyne detection, we obtain 



S{u, v) = 
with 



e 2^ 



1- z 



e i ^2;(l _ z) 



of the non-exponentially decreasing trajectories jumps 
to the ground state at exponentially distributed, finite 
times. For heterodyne or homodyne detection, however, 
it is not a priori clear whether the conditioned two-level 
system may become fully excited or deexcited at finite 
times as a result of the measurements. 

While deferring a rigorous verification to Sec. |Vl we 
will now argue heuristically for the attainability of the 
extremal values Ct = and Ct = 1, and in the process 
obtain an analytical expression for the excitation proba- 
bility Pt„<oo(2/)- 



1. Attainability of zero and unit excitation 

We suppose that the upper level population Ct = e << 
1 at time t. To attain the value in the next infinitesimal 
time step, Ct+dt = 0, the measurement must induce a 
change in Ct through a Wiener increment of size dWt ~ 
— r according to Eqs. (HI) and ([6]), where fc = for 



heterodyne detection and k ^ 2 for optimal homodyne 
detection. In the limit e J, 0, where the excited state 
population becomes infinitesimally close to 0, this would 
require an arbitrarily large Wiener increment, dWt — >■ 
oo. Consequently, the probability that the excited state 
population of an atom with a non-zero initial excitation 
y vanishes after a finite time is zero, and the first hitting 
time ro(j/) must be infinite with certainty. 



PTo=ooiy) = 1 



for y e ]0, 1] 



(16) 



For completeness, we note that a similar argument does 
not apply for the attainability of the excited state, since 
the diffusive term vanishes more slowly when Ct ap- 
proaches unity. If Ct = 1 — e is close to 1, the Wiener 
increment required for transitioning to Ct+dt = 1 after dt 
is dWt — >■ + 7"^), which may occur with finite 

probability in the limit e J, 0. 

Even though Ct does not reach in finite time, the 
emitter decays with certainty to the ground state in the 
long time limit. 



lim C, = 



(17) 



S{0,b)<oo, 5(a, l)<oo for a, & G [0, 1], (15) 

where erf(z) = e^^ dz' is the error function with 

erf(O) — and erf(oo) = 1. 



B. Attainability of ground and excited states and 
excitation probability 

While the unobserved, exponentially decaying emitter 
never exceeds its initial excitation and only reaches the 
ground state asymptotically, i.e. after infinite time, dur- 
ing continuous photon count detection, a fraction Cq = y 



2. Excitation probability 

The above analysis yields an explicit formula for the 
probability that Ct attains any level population u G [0,1] 
in finite time, 

r|i^ forO<y<u<l 

[l forO<it<?/<l 

where the first case follows directly from Eq. (fT2|) by 
putting b = u and taking the limit a \, Q upon using 
Eq. ([H]). We note from Eqs. (HH) and that this 
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limit and the limit u '[ I are well defined, and the for- 
mula also holds for u = 1. The second case follows from 
the fact that Ct tends continuously towards zero and con- 
sequently passes all points in ]0,y] with unit probability. 

For optimal homodyne detection, S{0, 1) < oo accord- 
ing to Eq. (jl5l) . and there is a non-vanishing probability 
that the upper level population Ct of the emitter (ini- 
tially prepared in any superposition different from \g)) 
attains any value u in ]0, 1] in finite time. This is il- 
lustrated in Fig. 2] for different initial level populations 
Co = y. In particular, the emitter has a finite proba- 
bility to become fully excited as a result of the optimal, 
continuous homodyne measurements. 

For heterodyne detection, S{0, 1) = oo from Eq. p4|) 
and, even though Ct may attain any value u G ]0, 1[, the 
excitation probability PT„<oo(y) drops to zero as u f 1, 
and the emitter has vanishing probability to become 
fully excited as seen in Fig. |31 More specifically, since 
Ei{^^) = -\n{^)-j+Oil-u), where 7 = 0.5772... 
is the Euler-Mascheroni constant, we see from Eq. ([18]) 
that Pt„ <oo(2/) varies as — i f^^^^^-^ in the limit w f 1- 
For not too small y, examination of this asymptotic ex- 
pansion (see Fig. m shows that PT^<oo{y) is finite for a 
wide range of values of u, while dropping steeply to for 
u very close to unit excitation. 

Finally, we note from Fig. |4] that compared to the opti- 
mal homodyne detection case, an atom subject to hetero- 
dyne monitoring is generally much less likely to become 
excited. 




FIG. 4. The probability Pt„<oo(j/), given in Eq. ([TSJ, tliat 
the upper level population Ct attains a specific value u above 
its initial value y = Co = 0.3, 0.5, 0.7, 0.9 (from left to right) 
for heterodyne (solid curves) and optimal homodyne (dashed 
curves) detection. The y- values are the same as in Fig. [2] The 
smaller figure presents a magnified view of the region near unit 
excitation probability. Ct only attains the value 1 with non- 
vanishing probability for optimal homodyne detection while 
Ct never reaches 1 for heterodyne detection. 



V. CHARACTERISTIC TIMES FOR 
HETERODYNE AND OPTIMAL HOMODYNE 
DETECTION 

A. First exit time and rigorous justification of 
attainability 

We will now further characterize the dynamics of the 
observed emitter by calculating the average of the first 
time Tafi{y) that the conditioned excited state popula- 
tion Ct, commencing from Co = y € [a, hits either a 
or 6. As a by-product we obtain a rigorous justification 
of the results on the attainability of zero and unit excited 
state population in Sec. IIVI 

In order to obtain a deterministic expression for 
E{Ta,biy)) = Ta,b{y), we abbreviate Ta,b{y) = T{y) and 
observe that the probability that the trajectory Ct, start- 
ing at time t — Q from y S ]a, 6[, will reach either a or 6 
after an infinitesimal time dt is vanishing. Thus by the 
Markov property of Ct, the first mean exit time of ]a, b[ 
at time dt is, conditioned of the value of C^t, given by 
T(Cif), and it is clear that E(r(Cdt)) = T{y)-dt, where 
the mean E(. . .) is taken over all possible realizations of 
Celt commencing from Co = y. 

Using the Ito rules as in Eq. we find 

T(y) =E(T(Crft)) +dt 

= Tiy) + T\v)A{y)dt + ]^t" {y)B{yfdt + dt, 

and we conclude that Ta.b(y) solves the differential equa- 
tion 



LyTa^iy) = -1 



(19) 



subject to the boundary conditions Ta,6(a) = Ta,bib) — 
0, and with the differential operator Ly given by the same 
expression Eq. pT|) as in the equation for Pr^KTAy)- See 
Ch. 15.3 of 13] and Ch. 5.2.7 of 14!] for further details 
on the derivation. 

We can verify that the solution of Eq. (|19p is 



Ta,b{y) 



S{y, b)R{a, y) - S{a, y)R{y, b) 



S{a,b) 

with y G [a, b], and < a < < 1, where 
R{u,v)= I Q{z)K{z)dz, K{z)^ 



Q{z')B{z'Y 



dz' 



and S and Q are specified in Eq. p^ . 



Heterodyne detection. 



For heterodyne detection, we obtain 



1 1 

K{z) = -e" 
7 

R{u, v) — — 
7 



2 1 
2~- + - 



In 



Z Z' 

- 2z 



(20) 
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with limiting behavior 

i?(0,6)=oo, i?(a,l) = oo for a, 6 e [0,1] (21) 

2. Optimal homodyne detection. 
For optimal homodyne measurements we obtain 
— (2z-l)e^ + ^erfi(^- 



i^(^) 



2z 
(22) 



u) = — rin(z) 

7 



2z 




1 - z 



-erfi 



2z 



where erfi(z) = dz' is the imaginary error func- 

tion with erfi(O) = and erfi(oo) — oo. We have not been 
able to evaluate the last integral in R in terms of stan- 
dard functions, but we observe the limiting behavior as 
the integrand in R remains bounded on [0, 1] 

i?(0,&) = cx), i?(a, l)<oo for a,& e ]0,1], (23) 



3. Attainability of zero and unit excited state population 

We are now in the position to provide stringent verifi- 
cation of the results in Sec. IIVI According to Thm. 1 in 
Ch. 5.21 of [H and Lemma 6.2 in Ch. 15.6 of [13], we 
have 

Theorem 1 

IfPn<To{y) > (i.e. ifS{0,b) < oc), thenToMy) = 
oo^R{0,y) = oo ^ PTo=oo{y) = 1- 

Theorem 2 

If Pn<T„{y) > and To^biy) = oo (i.e. if S{0,b) < 
oo and R{0,y) = oo), then either (i) To.b{y) — oo and 



Ct — or (ii) To^biy) < oo and Cj 



limt 



We verify from Eqs. (|T4l), (|T5|), ([^1), and (|23l) that the 
assumptions in Thms. [T] and [2] are fulfilled and Thm. [1] 
gives Eq. (IT51) . Eq. p7|) follows immediately if i) of 
Thm. [5] applies, while if ii) is true and the trajectory 
Ct hits b, since Ct is Markovian, the theorem may be 
applied again now with Ct proceeding from b, i.e. y ^ b 
and 6 — >■ &i > &. This argument may be repeated, and 
as we see from the first case of Eq. , the probability 
that Ct starting from bi hits bi+i > bi is strictly smaller 
than 1 and eventually i) of Thm. [5] will apply. 



B. Average excitation time 

Unlike the upper level population of the unobserved 
emitter, the conditioned trajectory may spend a finite 



amount of time above its initial excitation. More gen- 
erally, Ct, commencing from y G [a, 1], may spend one 
or more periods of time in an interval [i, r] C [a, 1] 
before exiting through a for the first time. We de- 
note this time T{y) and notice that its ensemble average 
E(r(j/)) = T{y) may be expressed as 



Taiv) 



f{Xt)dt 



(24) 



where f{z) = 1 for z e [t,r\ and f{z) — otherwise. 
Here Ta{y) is the first time Ct with Cq = y hits the level 
population a, as introduced in Sec. IIVI By ultimately 
taking the limit a | in Eq. ([24|) and employing that 
To{y) = oo, we will obtain an expression for the mean 
time spent by Ct in any interval [£, r] C [0, 1]. 

To obtain a differential equation for T(y), we 
break the integral up into two contributions, T(y) — 

E{J^' f{Ct)dt) E(/^^°^^^ f{Ct)dt). The first contribu- 
tion is just f{y)dt. And observing that the probability 
that Ct, commencing at time t = from y s]a, 1], reaches 
a during the infinitesimal time step dt is vanishing, we 
infer from the Markov property of Ct that the second 
contribution is equal to F,{T{Cdt)) , where the average is 
taken over all possible realizations of Cdt starting from 
Cq = y. Employing the Ito rules as in Eq. (j9|), we find 
that 



T{y)=f{y)dt + E{T{Cdt)) 



= f{y)dt + T{y) + T\y)A{y)dt + ^r{y)B{yfdt, 



and T(y) solves 



LyT{y) - -f{y) 



(25) 



subject to the boundary conditions T{a) ~ and 
■§ljT{y)\y=i = 0, and where Ly is given in Eq. (fTTj) . Since 
Ta{a) = 0, the first boundary condition is clear, while 
the second one follows by noticing that Ct is reflected at 
Ct — 1, and Ta{y) and thus T{y) should have a maximum 
at y = 1. For more details, see Ch. 5.2.7 of [13] for the 
special case f{y) = 1, and Prob. 15 in Ch. 15 of [l^] for 
the general case. 

With the prescribed boundary conditions, the solution 
of Eq. (US]) is 



Tiy)=Ria,y)-Sia,y)K{l), 



(26) 



where R and K are given by 



R{u,v) = 1^ Q(z)k{z)dz, k{z) = ^ Q(J)g(3)2 ^^'' 

and S and Q are given in Eq. (jl3p . The explicit formulas 
are lengthy, and we shall refrain from writing them here. 

By taking the limit a J, in Eq. (|26p and remembering 
that To(t/) = oo, we obtain an expression for the mean 
time the trajectory Ct resides in the interval [(,, r] C [0, 1]. 
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Since hits all level populations in ]0,1[ with non- 
vanishing probability and tends asymptotically towards 
0, this time must be infinite if £ = and finite otherwise, 
as an analysis of Eq. ([25)) also reveals. 

Putting £ = y and r = 1 in the limit a J, 0, T{y) yields 
the average time spent by Ct above its initial value Cq = 
y. We note that R{0,y) = and K{1) = K{1) - K{y), 
where K{\) is finite for both heterodyne and homodyne 
detection according to Eqs. ([20)1 and (l22t . and we find 

T{y) ^ S{0,y){K{y) ~ K{1)) (27) 

for y £ [0, 1] with r(0) = (from the boundary condi- 
tion) and r(l) = 0, since, in case of heterodyne detection, 
T{y) goes like — ^(1 — y)^ ln(l — y) when y t 1 according 

to the previously stated expansion of Ei{^-^). 

In Fig. [51 we plot the mean excitation time for both 
heterodyne and optimal homodyne detection, and ob- 
serve that it is notably larger for homodyne detection 
than for heterodyne detection for all initial upper state 
populations y e]0, 1[. 
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FIG. 5. The mean excitation time for heterodyne (solid 
curve) and optimal homodyne (dashed curve) detection given 
in Eq. ([?7)). 



C. Average time of decay 

In the Markov approximation assumed, monitoring 
does not modify the average decay properties of the emit- 
ter, it merely unravels the master equation solution into 
different quantum trajectories depending on the detec- 
tion scheme. Thus, on average Ct equals the uncondi- 
tioned upper level occupation probability Pee(i) given in 
Eq. ([2). Starting from an initial excitation Pee(O) — y, 
the exponential decay of pee{t) reaches an excited state 
population a at the time texpici) = 7~^ln^. Due to 
the fluctuations of individual trajectories, the first time 
the conditioned trajectory Ct hits a certain level a be- 
low its initial value Cq = y will be a fluctuating variable 



Ta{y), which may take both larger and smaller values 
than texp{a), as illustrated in Fig. [Tj 

By taking f{y) = 1 in Eq. (|25p . we obtain a differential 
equation for its average value F,(Ta{y)) — Ta{y) 

LyTaiy) - -1 (28) 

with the boundary conditions Ta{a) = and 
J-Taiy)\y^i^O. 

The solution of Eq. ^ is 

Taiy) = R{a,y)-Sia,y)K{l) (29) 

for < a < y < 1. In the limit a J, we get from Eqs. ([2T|) 
and ([23[) that To{y) = oo, as expected since the first 
hitting time of is infinite by Eq. p^ . For heterodyne 
detection, the expansion of i?i(i^) implies that both 

R{a, y) and S{a, y)K{l) diverge like — ^ \n{l — y) as y f 1; 
thus leaving Ta(l) finite. 

In Fig. [6l we compare the average exponential decay 
in Eq. ^ with the average first passage of excitation 
levels below the initial population for both heterodyne 
and optimal homodyne detection given in Eq. (12^1) . The 
figure shows that on average Ct conditioned on hetero- 
dyne detection reaches lower values faster than pee{t)- It 
is interesting to observe that homodyne detection, which 
leads to excursions to higher degrees of excitation, also 
shows the faster average approach to low excitation levels 
as compared to heterodyne measurements. 
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FIG. 6. Comparison of exponential decay (dotted curve), 
(t'Y, pee{t)) , with the average time to reach excited state 
probabilities a below y, {"yTaiy), a) , during heterodyne (solid 
curve) and optimal homodyne (dashed curve) detection. The 
initial excitation is pee{0) = y ~ 0.5. The two-level system 
conditioned on continuous, heterodyne or optimal homodyne 
measurements reach lower excitation faster on average than 
the unconditioned master equation solution. 
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VI. CONCLUSION 

The current interest in investigating the properties and 
achievements of conditioned dynamics are stimulated by 
the prospects of combining measurement back-action and 
feedback for control of quantum systems. We have in- 
vestigated the stochastic evolution of a two-level emitter 
subject to continuous homodyne and heterodyne detec- 
tion. In particular, we have found that the emitter may 
reach any excited level population, including unit exci- 
tation, as a result of in-phase, optimal homodyne mea- 
surements. Homodyne detection 7r/2 out of phase with 
the emitting dipole, on the other hand, cannot increase 
the excitation beyond its initial value, while homodyne 
detection with other local oscillator phases as well as het- 
erodyne detection may lead to any excitation probability 
strictly below unity. 

These results do not violate energy conservation since 
the initial state is assumed not to be an energy eigenstate. 
The results obtained are characteristics of the measure- 
ment process and of its ability to distinguish and unravel 



the dynamics into different conditioned trajectories. A 
projective measurement in the energy eigenstate basis 
would, indeed, yield the fully excited state with prob- 
ability equal to the initial excited state population. It 
is interesting, however, that heterodyne and homodyne 
measurements are based on the detection of the radiation 
emitted by the source, and that, in order to provide this 
field, the system has to loose energy. In a photon count- 
ing experiment, the detection of an emitted quantum of 
energy would be associated with a jump of the emitter 
into its ground state, while intervals of no counts would 
cause a continuous reduction of the excited state popula- 
tion [sj, and the excitation would not increase above its 
initial level. 

Conditioned dynamics has been demonstrated to be a 
useful and efficient way to prepare quantum states which 
are difficult or even impossible to reach by unitary time 
evolution, and our analysis of how even a very simple sys- 
tem can be prepared in an excited state by monitoring its 
decay, points to rich opportunities for further exploration 
of the achievements of different detection schemes. 
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